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We consider an Ising model where longitudinal components of every pair of spins have antiferro- 
magnetic interaction of the same magnitude. When subjected to a transverse magnetic field at zero 
temperature, the system undergoes a phase transition of second order to an ordered phase and if 
the temperature is now increased, there is another phase transition to disordered phase. We provide 
derivation of these features by perturbative treatment up to the second order and argue that the 
results are non-trivial and not derivable from the known results about related models. 
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I. INTRODUCTION 

Ising model with long range interaction is of interest in 
statistical physics mainly in connection with the entan- 
glement properties of the model and the quantum phase 
transitions displayed by it. A widely studied model of 
this type is the Lipkin-Meshkov-Glick (LMG) model ([l|, 
and references therein) which deals with a system of 
spin-i particles, each interacting with all other by the 
same strength and subjected to an external field. The 
general form of LMG Hamiltonian is 
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where S a (for a — x, y, z) are the components of total 
magnetic moment, 
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(2) 



N is the number of spins, J x , J y are the strength of 
interaction along X and Y directions for every pair of 
spins and T is the strength of an external field in the Z 
direction (s z = ±1/2). This model has been attacked 
in various ways (0] and references therein). One way is 
to ignore the non-commutation of the operators S a and 
treat them semi-classically : 

S x =Ssm6 cos ((>, S v = S sin(9 sin S z = <S cos 6 (3) 

where S is the total spin. Another way [3j is to use 
Gaudin Lie algebra and derive an exact solution in the 
form of Bcthe-like equations. 

In this paper, we consider a variant of LMG model 
which has been introduced recently and studied using 
quantum Monte Carlo method by Chandra, Inoue and 
Chakrabarti Q. One has an antiferromagnetic interac- 
tion of constant strength only in the longitudinal direc- 
tion, and an external transverse field. The Hamiltonian 
is 
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with J > 0. The LMG Hamiltonian reduces to this when 
either of J x and J y is zero and the other is negative. This 
Hamiltonian shows interesting phase transition proper- 
ties as a function of field strength T and temperature. 
Thus, when T is zero, the system is in a highly degen- 
erate state with zero longitudinal (total) magnetic mo- 
ment at zero as well as non-zero temperature. But at 
zero temperature, as soon as a small T is switched on, 
the system gets oriented completely in the transverse di- 
rection. When temperature is increased, keeping T non- 
zero, the system undergoes a transition from an ordered 
to disordered state at a certain critical temperature. The 
objective of this paper is to study the phase transition 
properties of the Hamiltonian TL using perturbative treat- 
ment. The first order term is zero and the second order 
term is calculated exactly. In the next section we shall 
derive the eigenvalues of Ti. for small values of T and in 
the last section present some features of the result. 

Before we conclude this section, we point out that the 
solution for the Hamiltonian H does not follow from the 
known treatments of LMG Hamiltonian and hence the re- 
sults of this paper are non-trivial. Since the commutator 
of the two terms in TL, 
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has expectation value ~ N when the spins are oriented 
in the transverse direction, we cannot make the semi- 
classical approximation described in Eq. (3). On the 
other hand, the parameter region covered here, although 
excluded by some discussions on LMG model Q , are in- 
deed included in the exact solution of Ortiz et. al. p|. 
However, the solution provided there is not in the form 
of a closed expression. We believe therefore, that the 
solution presented here are non-trivial. 



II. PERTURBATIVE TREATMENT OF H 

We shall study the Hamiltonian 7i for small values of 
r, by treating 
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as the perturbation over the (unperturbed) Hamiltonian and treat 



The eigenvalues of Ho are obviously 
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with degeneracy 
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where, 
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Let us call the set of eigenstates with eigenvalue E^ 
as M.. Since H p operating on any spin distribution 
\ a) e M gives a state with a different S z , the quan- 
tities (f3 \ Hp | a) is zero for all | a) , | (3) E M . The first 

order perturbation correction E$ is hence nil and the 
dominant correction comes from the second order correc- 
tion Ein ' . 

(2) 

The second order perturbation correction Em to the 
eigenvalues are the eigenvalues of a matrix P whose 
elements are given by 
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where a), | /?) but | I) M. We shall now see 

that the eigenproblem of P cn be solved exactly First 
consider the diagonal element P aa . The state | a) has 
(N/2)+m up-spins and (N/2) — m down-spins. When the 
operator H p flips an up/down spin E^ becomes (to + 
1) J, so that 



p = 



NT 2 4to 2 + TV 
4 J 4m 2 - 1 



(12) 



As regards the off-diagonal elements P a f3, we note that 
it is non-zero when and only when the spin-distributions 
a) and | j3) differ in precisely two unlike spins (since 
these two states must have the same S z ). Thus, the two 
spin-states will be like : 



\a) :\ 
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Then | 1} can be either | h 1 ) (with E\ ' = 

(m + 1) J) or | ) (with £, (0) = (m - 1) J). 

Hence, -P Q/ 3, if non-zero will have the value 
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To proceed, let us now consider the Hamiltonian 

J_ {sz) 2 h 
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, h r 
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as perturbation on the (unperturbed) Hamiltonian 
jj (S z f which is nothing but the H of Eq. (5). The 
first order perturbation matrix P' has diagonal elements 

PL = (a\HL\a) 
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The off-diagonal elements P' a g will be h/N if | a) and 
| /3) differ in precisely two unlike spins and will be zero 
otherwise. Thus, we get the correspondence 
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where 1 is the VL m x fi m unit matrix. 

It is trivial to diagonalise the Hamiltonian Ti' by 
rewriting it as, 



v = --LJl (s'-f 



(18) 



(where S is the total spin operator) and noting that the 
eigenstates of this Hamiltonian is labelled by to and the 
total spin S. Using the usual relationships 

S 2 | S, to) =S(S + 1) I S, to) 



N 

5 = 0,1,2,. 

and 



m = -<S,-5 + l,---5-l,5 
we observe that the eigenvalues of P' are 
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From Eq.(17) we can finally conclude that the unper- 
turbed eigenstate of Ho with eigenvalue (J/N)m 2 splits 
up under second-order perturbation by H p to eigenstates 
with eigenvalues 

NT 2 
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The degeneracy of this state can be easily seen to be 
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III. DISCUSSION 



state is, 



Having solved the eigenproblem of TL upto second order 
in r, we can now observe the following features : (i) At 
zero transverse field, the system has Hamiltonian 7io and 
is in a state of m = for zero or non-zero temperature, 
(ii) At zero temperature, when a small (transverse) field 
is turned on, the ground state energy (correct upto V 2 ) 
becomes, 
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with S = N/2. The susceptibility per spin, 
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diverges in the thermodynamic limit, indicating a sec- 
ond order phase transition at T = 0. (iii) At non-zero 
temperature T, the free energy (correct upto T 2 ) is, 
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where k B is Boltzmann constant. The value of S which 
minimises this expression is hence the solution of the 
equation 



x = tanh(ax) 



(25) 



where x = 2S/N, and a = N 2 T 2 / (AJk B T). To arrive at 
this relation we have assumed N, S and (N/2 — S) to be 
large and used Stirling's approximation for the factorials 
of Eq.(21). This shows that there is a mean-field type 
phase transition at a critical temperature 
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Two aspects of this phase transition are to be noted, (i) 
Since the expectation value of (S x ) in the unperturbed 



at T > T c 
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there is disordered phase with 
5 = and ((S X Y) = 0, while for T < T c , we have an 
ordered phase with S ^ and ((S x ) 2 } / 0. (ii) The ap- 
pearance of N in the expression of T c is uncomfortable 
at the first sight. Note that for the perturbative treat- 
ment to be valid, the perturbation correction must be 
less than the spacing between two adjacent eigenvalues 
of the unperturbed Hamiltonian. According to Eq. (4) 
this needs, 



r< J/N 



so that 
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Thus, all we can claim is that when the condition (27) is 
fulfilled there will be a temperature-induced phase transi- 
tion at a temperature given by Eq. (26). In other words, 
if we rewrite the Hamiltonian of Eq. (4) as 
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(where b = TN/J), then for 6 = 0, there is no critical 
temperature but when b has a small non-zero value, there 
is a critical temperature given by T c = b 2 /A and this 
conclusion is true for all values of N. Our treatment 
cannot predict anything when b is not small. In the finite 
size study of Chandra, Inoue and Chakrabarti [3| the 
value of b is 100 and hence the fact that this study does 
not observe a thermal phase transition does not present 
a contradiction with our results. 
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